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In this paper we investigate the behaviour of the specific heat around the critical point of the Ising 
model in dimension 5 to 7. We find a specific heat discontinuity, like that for the mean field Ising 
model, and provide estimates for the left and right hand limits of the specific heat at the critical 
point. We also estimate the singular exponents, describing how the specific heat approaches those 
limits. Additionally, we make a smaller scale investigation of the same properties in dimension 6 
and 7, and provide strongly improved estimates for the critical termperature Kc in d = 5, 6, 7 which 
bring the best MC-estimate closer to those obtained by long high temperature series expanions. 
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INTRODUCTION 

The Ising model in dimension d = 5, being strictly 
larger than the upper critical dimension dc = 4 of the 
model, has been studied by many authors, and has been 
the focus of a long running debate regardi^ its finite size 
scaling behaviour. We refer the reader to [l| for a discus¬ 
sion of that topic. However, even regarding the infinite 
size limit of the model there are still interesting open 
questions. For d > dc it is rigorously known that on the 
d-dimensional hypercubic lattice the critical exponents of 
model takes their mean field values. This was first proven 
in [2-01 and more recently the newly developed rigorous 
lace expansion for the Ising model [6| made it possible to 
give a unified proof of these results by a single method 
[7| . For the specific heat the critical exponent a = 0 and 
the results of [2] also show a stronger result, namely that 
the specific heat is bounded at the critical point. How¬ 
ever, the values of the critical exponents are not the only 
properties which are characteristic for the phase transi¬ 
tion in the mean field version of the model, there is also 
a discontinuity in the value of the specific heat at the 
critical point. In fact, historically this discontinuity was 
once noted as one of the first signs showing that mean 
field theory does not give a correct description of phase 
transitions in dimension 2 and 3. Hence it is natural to 
ask if there is a similar discontinuity in the specific heat 
for the Ising model above the upper critical dimension. 

This question is also natural from another point of 
view. One way of realizing the mean field version of the 
Ising model is to view it as the thermodynamic limit of 
the Ising model on finite complete graphs. The model 
on complete graphs has been studied rigorously in great 
detail by mathematicians, both in the usual Ising form 
and the equivalent Fortuin-Kateleyn random cluster rep- 
resentation, [2, Q . Numerically it has been observed [i2| 
that for d = 5 the model on finite lattices with periodic 
boundary conditions displays the same scaling behaviour 
inside the critical scaling window as the Ising model on 
a complete graph. Hence it is also natural to ask if the 
behaviour of the thermodynamic limit of the model on 


the complete graph and the hypercubic lattice will also 
show the same type of behaviour at the critical point, in 
particular if both have a discontinuous specific heat and 
how large the jump at that discontinuity is. 

In order to study this question we have done Monte 
Carlo simulation of the Ising model on hypercubic lat¬ 
tices with periodic boundary conditions, with the main 
effort for d = 5 but with some data for d = 6, 7 as well. 
Using these data we first give improved estimates for the 
critical temperatures in these dimensions. Next we find 
that there is a jump in the specific heat and give esti¬ 
mates for the left and right hand limits of the specific 
heat at the critical temperature Kc- We also estimate 
the singular critical exponents, describing how the spe¬ 
cific heat approaches the limit values. As mentioned in 
0 the singular exponents, unlike the critical ones, are 
not expected to necessarily have the same value on the 
low and high-temperature sides of the critical point, and 
we find that their values are quite distinct. Finally we 
also note that as d increases the behaviour at the critical 
points seems to be approaching that of the mean field 
limit, as expected. 

The structure of the paper is as follows. After some 
definitions we first give a derivation of the specific heat 
for the complete graphs, and use it to give a description 
of the specific heat for the mean field limit which is more 
detailed than the usual one. Next we present our nu¬ 
merical data for d = 5, first for the critical temperature 
Kc and then for the specific heat near Kc- After that 
we give a brief description of the corresponding results 
for d = 6,7, and finally we give some discussion of the 
observed results. 


DEFINITIONS AND DETAILS 

For a given graph G on N vertices the Hamiltonian 
with interactions of unit strength along the edges is 
% = — J2ij SiSj where the sum is taken over the edges 
ij. As usual the coupling K = l/ksT is the dimen¬ 
sionless inverse temperature and we denote the thermal 
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equilibrium mean by (•••). 

The we call the critical coupling itTc, and denote its 
normalised form by £ = {K — K^jK^ and the rescaled 
version hi = \fN{K — Kc)lKc. As usual the magneti¬ 
sation is M = Si (summing over the vertices i) and 
the energy is E = SiSj (summing over the edges ij). 
We let m = M/N, U = E/N and U = (U). The specific 
heat is defined as 


—d'^ logZ 
dTdK N 


“/V 


(( a 2 ) - (Ef) 
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When the underlying graph is a d-dimensional grid 
graph of linear order L with periodic boundary condi¬ 
tions we mean it simply to be the cartesian product of 
d cycles on L vertices, so that N = L‘^. When we re¬ 
fer to the complete graph we mean the graph where all 
pairs of vertices are connected by an edge, thus having 
(^) “ N(N — \)I2 edges. For d = 5 we have collected 
data using Wolff-cluster updating for L = 16, 20, 24, 32, 
40, 48, 56 and 64. The number of measurements at each 
temperature near Kc ranges from ca 15000 for L = 64 
to more than 100000 for L = 16, 20, 24,32. We will also 
re-use some extremely detailed data from [l^ for L = 6, 
8, 10 and 12. 


THE ISING MODEL ON THE COMPLETE 
GRAPH 

Recall that the complete graph Gn on N vertices is 
the graph with N vertices in which every pair of of dis¬ 
tinct vertices is joined by an edge. The limit as A —>■ oo 
of the Ising model on Gn corresponds to the usual mean 
field Ising model. In order to be able to make a detailed 
comparison with the d-dimensional Ising model we will 
now derive an expression for the specific heat of the mean 
field model in a neighbourhood of the critical point, in¬ 
stead of the more common textbook version which only 
gives the jump exactly at the critical point. 

First note that it is an exercise to show that (E) = 
(1/2) — N/2 and, more importantly, 


Pr(M = —2) = Pr(M = 0) = Pr(M = +2) (for even 
N) and thus 


Kc 
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constitutes an effective Kc- The appendix of Ref. [l^ 
provides detailed information on the moments of this 
magnetisation distribution and we will apply this to get 
information on the energy moments. We begin with the 
case oi q = N/{N -I- 2) — 2 «;/iV^A, With q = ex.p{—2K) 
this corresponds to k = \/N{K — Kc)/Kc + 0{1/N), 
i.e. we move around inside the scaling window with the 
temperature parameter k. Using Lemma A6 of Ref. [l^ 
(after setting a = —2k) we can now easily obtain the 
asymptotic form of the £th moment as 


OO 

N3e/^2^ f |a;|^i?(K,x)da; 

{\Mf) -- (5) 

f R{K,x)dx 

— OO 


where R{k,x) = exp(2Kx^ — 4x'^/3). Plugging this into 
Equation m and evaluating the integrals we can obtain 
a formula for C. However, in the special case k = 0, we 
get the very simple 

= i - 

The local maximum of C can now be computed numeri¬ 
cally to lie at K* = 2.2568473919660... and the value at 
this point is Cmax = 1.6572974585496.... We note also 
the limits limK->.oo C(k) = 3/2 and limK->.-oo C(«^) = 0. 

To continue with the case outside the scaling window 
we set q = {n — 2e)/(n -|- 2) which, since q = exp(—2Ar), 
gives us e = {K — Kc)/Kc + 0{\/N), our normalised 
temperature. In the high-temperature case, i.e. for e < 
0, we get from Lemma A9 (setting a = —2e) of Ref. [I^ 
that 

OO 

NiG2‘^ f |x|^exp(2£x2)dx 

{\Mf) = -- (7) 

J exp(2ex^)dx 


var(F;) = ivar(M2) = 1 ((M^) - (2) 

This will come in handy when we compute C{K) = 
\inlN-^ooC{K, N) where C{K,N) = K^vai{E) /N. 

It is shown in Ref. [1^ that the magnetisation distri¬ 
bution at coupling K for a complete graph is 

Pr(M = A - 2fc) = , 0<fc<A (3) 

where q = exp(—2Ar). Since = N — 2k) = I 

this implicitly defines i)). When q = N/{N + 2) the 
distribution is precisely flat in the middle, i.e. with 


so that C{e,N) = 0{1/N) and thus C(e,oo) = 0. This 
could be interpreted as 0“ = I. The case £ = 0 was 
treated above as k = 0. 

The low-temperature case £ > 0 is a little more tricky. 
Let n = (|to|), where 0 < p, < I, denote the normalised 
(spontaneous) magnetisation and note that the magneti¬ 
sation distribution has a peak at Mpeak = ifJ-N having 
width 0{\fN). Moving x\fN magnetisation steps away 
from Mpeak we are at the new magnetisation where 
\Mx\ = 2^fN\x + p.'/N/2\. Lemma AI3 says that the 
ratio Pr(Ma;)/Pr(Mpeak) is asymptotically 

R{^i, x) = exp |2x^ (^ ^2 1 ^ I (8) 
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FIG. 1: (Colour on-line) C{£,N) versus e, with e = {K — 
Ka)/Kc for finite N = 32, 64, 128 256, 512, 1024 (increasing 
blue curves) and the complete graph limit case N ^ oo (red 
thick curve). The red dot is the limit at e = 0. 


and the £th moment then becomes 


(|M0 


iV^/22^ / 




R{n,x)dx 


f R(/i,x)dx 


(9) 


Using Equation ([2]) the specific heat limit, expressed 
in /i, collapses into the simple form 


C(m) = 




^ — 1) atanh(^) 


( 10 ) 


Next, Lemma All [13|, after setting a = —2e, says that 
/i depends on e as asymptotically 


In combination with Equation (US this defines implicitly 
the limit specific heat in terms of e. Taking the compo¬ 
sition of the series expansion of Equation (US and the 
inverse series expansion of Equation (HU we obtain at 
last 


C{e) 


3 12 438 n 432 ., 

- £ H-- £^ + 

2 5 175 175 


166104 

67375 


£^ + --- ( 12 ) 


Plotting the numerical evaluation of m and (HU we 
get the Figure [T] where the limit and some finite cases 
are shown. 


THE 5-DIMENSIONAL CASE 

We now come to our Monte Carlo results for d = 5. As 
noted in [iU| some of the indicators used by other authors 
to study the critical behaviour for d = 5 are very sensitive 


to the exact value of the critical temperature K^- With 
this in mind we will first present a new way of obtaining 
highly precise estimates for Kc and use it to derive the 
estimate which we will use in our later analysis. 


An improved method for estimating Kc 


Our improved estimate of suitable for d > 4, is 
based on a careful study of the magnetisation distribu¬ 
tion, i.e. Pr(M). The approach is simple but assumes 
that all measurements of M were stored for each tem¬ 
perature during the sampling process. After normalis¬ 
ing these values as a; = we put them in bins 

of reasonable width, in our case 0.20, thus giving us a 
histogram. There are of course several different binning 
methods to choose from, but for simplicity we have cho¬ 
sen to use a fixed bin width which is roughly what the 
Freedman-Diaconis method (twice the interquartile range 
divided by the cube root of the number of measurements) 
prescribes when the distribution is near Kc{L) (see be¬ 
low) for the weakest data set (i.e. for L = 64). 

The simple distribution density function fix) = 
4>o exp((^ 2 a;^ + 4>ax'^) is then fitted to this histogram. 
Since 4)2 for all intents and purposes depends linearly 
on K inside the scaling window (see Fig. 3 of 0 ) , we 
fit a straight line to the data points (at least seven) on 
the interval corresponding to —0.7 < 4>2 < 0.7 and solve 
4>2iK) = 0. This point constitutes an effective criti¬ 
cal temperature Kc{L) scaling as Kc{L) — Kc cx. 
see 0. 

Ideally the density function fix) should also contain a 
correction factor (1 -|- X 2 X^ + X^x'^ -I- Aqo;® -I- ...) but the 
coefficients Xi will vanish with increasing L. For L > 16, 
especially near KdL)^ they will not contribute signifi¬ 
cantly to fix) and can in any case not be discerned with 
the data we rely on here. See 00 for a consider¬ 
ably more detailed study of the scaling behaviour of the 
magnetisation distribution. 

In Figured] we show KdL) versus 1/L^ together with 
an inset showing the M/N^/^-distribution for L = 32 at 
different values of K and another inset showing how 4>2 
depends on K for the different system sizes. A line fit 
gives that KdL) = 0.11391498(2) - 0.0654(2) L-^. The 
coefficients and their error estimates are here based on 
the median and interquartile range of the fitted coeffi¬ 
cients when deleting each data point in turn from the 
line fit. The estimate Kc = 0.11391498(2) is within the 
error bars of earlier estimates 0,0 but adds another 
digit to the accuracy. This technique for estimating Kc is 
quite robust to variations in the various parameters. For 
example, changing the distribution bin widths to 0.15 or 
using 4>2 data for —0.6 < 4>2 < 0.6 keeps the resulting Kc 
within the stated error bars. 
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FIG. 2: (Colour on-line) Kc{L) versus 1/L® for L — 16, 20, 24, 
32, 40, 48, 56 and 64. The fitted line is 0.11391498 - 0.0654a:. 
The lower inset shows the fitted parameter (j )2 versus K for the 
same L. The upper inset shows the normalised magnetisation 
distribution for L = 32 at if = 0.1139, 0.113915, 

0.11393 and 0.113945 resulting in 02 = -0.251, 0.036, 0.323 
and 0.621 respectively of the fitted fix) (red curves). 
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FIG. 3: (Colour on-line) CiK,L) for L = 4, 6, 8, 10, 12, 16, 
20, 24, 32, 40, 48, 56 and 64. 


Specific heat discontinuity 

Consider Figure [3] where we plot the specific heat for 
4 < L < 64 for a wide temperature range. Clearly there 
is an envelope curve containing the limit specific heat. 
From the individual CiK,L) functions we extract the 
limit function C{K,oo) = limi_).oo C(Ff, L) from points 
where the function for increasing L agree. We thus 
assume that there is a Kj^iniL) such that if L' > L 
and K > Knun{L) > Kc then CiK,L') = C{K,oo). 
Analogously we assume there is a K^^^iL) such that 
C{K,L') = C{K,oo) when 0 < AT < 
and L' > L. As an example, taking for example Li = 16 
and L 2 = 20 we note that C(Ar, 16) = C(Ar, 20) when 
K > Ffinin = 0.1148, where ATmin of course depends on 



FIG. 4: (Colour on-line) C{K,L) versus A for A < Ai„ax(T) 
for L = 6, 8, 10, 16, 20, 24, 32, 40, 48, 56 and 64. The 
curve consists of more than 500 points and error bars are not 
shown. The inset shows a zoomed-in version with 3rd order 
interpolations through all data points. Data for increasing L 
start to deviate from envelope cnrve as we move closer to Kc- 

the chosen Li and A 2 . On the high-temperature side 
of Kc we find that C(A, 16) = C(A, 20) when A < 
ATmaxClfi) = 0.1120. Thus we treat the measured data 
for C(A, L) as the asymptotic C(A, 00 ) when A > 16 and 
K < 0.1120 or A > 0.1148. An increasing sequence of 
pairs of Ai, L 2 gives a sequence of Amin and Amax that 
both approach Kc- The individual Amin and Amax were 
found by simply comparing pairwise plots of C{K, A). 

In Figure m we show the individual C(A, A) for A < 
Kniax{L) pieced together into one plot for a range of 
A and Figure [S] shows the coresponding data for A > 
Amin(A). Their insets shows the data without remov¬ 
ing the finite size behaviour and clearly demonstrate the 
presence of a limit enveloping curve. Having removed any 
finite size effects, such as the local maximum for each A, 
the remaining points are in effect estimates of the asymp¬ 
totic C(A, 00 ) for 0 < A < Amax(64) = 0.11388 and 
A > Ami„(64) = 0.11394. 

In order to estimate the left- and right-limit we take 
the limit curves and fit a simple expression of the form 

Aq -|- Aix^[1 B\x -|- B2X^ -|- ...) ( 13 ) 

to C(A, 00 ) where x = \K — Kc\/Kc- We will use ±- 
superscripts to denote the left- and right-limit as a; —>■ 0. 
This provides left- and right-limits (A(C and A({"), the 
dominating correction term exponents (0“ and 9 ^) and 
a sequence of correction terms. Following [HI we call 9 
and 0+ the singular exponents of the specific heat, since 
they describe the behaviour of the singular part of the 
specific heat. We are not aware of any prescribed form 
of the correction terms from earlier studies so these will 
simply be the effective terms. 

Using Mathematica’s built-in FindFit-function we fit 
the high-temperature limit curve to (1131) using both one, 
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FIG. 5: (Colour on-line) C{K,L) versus K for K > Kmin{L) 
for L = 6, 8, 10, 16, 20, 24, 32, 40, 48, 56 and 64.Tlie curve 
consists of more than 500 points and error bars are not shown. 
The inset shows a zoomed-in version of 3rd order interpola¬ 
tions through all data points. Data for increasing L deviate 
from envelope curve as we move closer to Kc- 

two and three correction terms and find excellent agree¬ 
ment in the resulting values of the singular exponent 9~, 
measuring 9~ = 0.40(1). The first two coefficients 
and also strongly agree when adding more correction 
terms. However, having first established a strong candi¬ 
date exponent we now simply fix this to 0~ = 0.4 and use 
(nsj, again trying one, two and three correction terms. 
Based on this we find an effective fit 

C-(x) = 0.1697(2) - 0.231(l)a;° '‘°(l - 0.26(l)x) (14) 

where x = —e = {Kc — K)/K and 0 < K < Kc- Adding 
more correction terms does not improve the fit. The error 
bars reflect how the coefficients change when adding one 
or two more terms. 

Repeating this exercise for the low-temperature side 
the FindFit-function suggests 0+ = 0.60(2) and again 
the leading coefficients agree using one, two and three 
correction terms. Setting 0+ = 0.60 gives us the effective 
fit 

C+(x) = 2.040(1) - 2.58(l)x° ®°(l - 0.36(2)x) (15) 

where x = e = {K — Kc)/K and K > Kc- As before, 
the error bars reflect how the coefficients change when 
adding correction terms. We now put m and m to 
the test by taking log-log plots of the measured C{K, oo) 
when subtracting the respective limit A^. 

Beginning with the high-temperature case, in Figure [S] 
we show log(0.1697 —C(Ar, oo)) versus log(x), where x = 
{Kc — K)/Kc, together with C~{x) of (ITTl) (black curve) 
and the asymptote 0.1697 — 0.231x° "‘° (red line). The 
rather small error bars suggest a good quality of the fit. 
Analogously, on the low-temperature side, we show in 
Figure [7] log(2.04 — C(Ar, oo)) versus log(x), where x = 


]n|A„-C| 



FIG. 6: (Golour on-line) log(0.1697 — C(A, oo)) versus log(x), 
with X = {Kc- K)/Kc, for L = 6, 8, 10, 16, 20, 24, 32, 40, 
48, 56 and 64 (larger L to the left in the figure) together with 
the fitted C~{x) (black curve, hard to see) of Eq. (I14II and the 
asymptote 0.1697 — 0.231x°''* (red line with slope 0.40). 


]n|A„-C| 



FIG. 7: (Golour on-line) log(2.04 — C(A, oo)) versus log(x), 
with X = {K - Kc)/Kc, for L = 6, 8, 10, 16, 20, 24, 32, 40, 
48, 56 and 64 (larger L to the left in the figure) together with 
the fitted (x) (black curve, hard to see) of Eq. (I15II and the 
asymptote 2.04 — 2.58x° ®° (red line with slope 0.60. 

{K — Kc)/Kc, together with C+(x) of (IT5)l (black curve) 
and the asymptote 2.04 — 2.58x°'®° (red line). For L = 
56, 64 the error bars are now quite pronounced. For the 
smaller L the error bars are considerably more benign. 

Finally we estimate the value of C{Kc,oo), which, of 
course, does not have to coincide with any of the A J. As 
it turns out this value is quite distinct from both limits. 
In Figure [5] we show a zoomed-in plot of C{K,L) over a 
range of L for three fixed AT-values, K = 0.11391498 (i.e. 
the estimated Kc), K = 0.1139148 and K = 0.1139152. 
The C{K, L) were found by interpolating the data points. 
As the plot demonstrates, there is a clear upwards trend 
in the values for K = 0.1139152 and a clear downwards 
trend for K = 0.1139148, whereas the middle value shows 





















6 




FIG. 8: (Colour on-line) C{K,L) versus for L = 16, 

20, 24, 32, 40, 48, 56 and 64 together with the fitted red line 
0.724 -I- 4.19a;, where x = . The points trending up¬ 

wards are for K = 0.1139152, the points trending downwards 
are for K = 0.1139148 and the middle set of points are for 
Kc = 0.11391498. The rectangle on the y-axis indicates the 
estimate C{Kc,oo) — 0.724(3). 

no clear trend. A fitted line on the points for L > 16 
suggests C{Kc, oo) = 0.724(3). The error bar of this value 
is obtained by allowing the value of K to vary inside the 
error bar of Kc (2 steps in the 8th digit) and repeat the 
line fit to the new points. 

The local maximum of C(A', L), see FigurejSl also takes 
its own limit value, i.e. Cmax = limL_>.oo niaxif C(Ar, L) 
does not coincide with the right-hand limit C~. In 
Figure IHl we show the estimated maximum for each L 
and the right-hand limit 2.04 found above (HSl). It ap¬ 
pears very unlikely that they should coincide for large 
L. The fitted line, based on L > 16, suggests a limit 
Cmax = 2.225(6) where the error bar is based on the 
variability of the constant term of fitted lines (with 
X = with one point removed from the data set 

L > 16. We estimate that the maximum is located at 
= Kc + 1.860(2)F^c/A®/^ i-e., at k = 1.860(3), 
with the error bar obtained as before by removing indi¬ 
vidual points for L > 16 when fitting a line through the 
origin (since we know ATmax —t Kc) ■ We will plot C versus 
K later. 

We can now make a comparison of the behaviour of 
the 5-dimensional model and that of the mean field case, 
as derived in the previous section. We first consider the 
scaling window, in Figure (TU] we show a plot of C{k) for 
a range of N and the limit case, together with our data 
for d = 5. As we can see that maximum specihc heat for 
the mean field limit is lower than the values for d = 5, 
but the general shape of the curves are nonetheless quite 
similar. 

Next we look at the thermodynamic limit. In FigurefTT] 
we show the specihc heat limit for both the complete 
graph and d = 5 in the same plot. As we just noted, the 
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FIG. 9: (Golour on-line) maxK C{K, L) versus , for L = 

16, 20, 24, 32, 40, 48, 56 and 64 together with the fitted red 
curve 2.225 — 48.6a;, where x = The rectangle on 

the y-axis indicates the limit estimate Cmax = 2.225(6). The 
point at y = 2.04 is the right-limit C'*’(0) = 2.04 of Eq. (1151) . 
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FIG. 10: (Golour on-line) C versus k, with k = '/N{K — 
Kc)IKc for finite N = 32, 64, 128 256, 512, 1024, 2048 
(increasing blue curves) and the complete graph limit case 
A —^ oo (red thick curve). The upper set of points shows the 
same for 5D sampled data points for L — 16, 20, 24, 32, 40, 
48, 56 and 64. Error bars are shown but of the same size as 
the points. 

value for the mean field are lower than those for d = 5 
when we are sufhciently close to e = 0. We can also see 
the difference in the singular exponents between the two 
models, with the mean field case approaching the line 
e = 0 at an angle and the d = 5 case instead approaching 
it tangentially. 

THE MODEL IN DIMENSIONS 6 AND 7 

As the dimension d increases we should see the specific 
heat approach that of the complete graph. We also col- 
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FIG. 11: (Colour on-line) The limit specific heat C{e with 
e = {K — Kc)IKc for the complete graph (red, y = 1.5, y = 
0.407 and y = 0 at y-axis) and the 5D case (black, y = 2.04, 
y = 0.72 and j/ = 0.17 at y-axis). 

lected data for d = 6 and d = 7 and tried to estimate the 
singular exponents 6~^ and 6~. However, these data rely 
on considerably smaller systems; L = 4,6, 8, 10,12,16, 20 
for d = 6 and only L = 4, 6,8, 10 for d = 7. The process 
is the same as we used above for d = 5 and we will simply 
state the resulting estimates of the various parameters. 

For d = 6 we estimate Kc = 0.0922982(3) and 
C{Kc,oo) = 0.58(1). This value of Kc deviates some¬ 
what from the older Monte Carlo estimates, as surveyed 
in [l^, which have tended to be close to 0.09229, but 
agrees well with the more recent estimate 0.092298(1) 
[l9|, coming from the longest series expansion results to 
date. The limit specific heat is 

{ 1.833- 2.61e°-^5('i_o.57£:) e>o 

0.58 e = 0 (16) 

0.0927 - 0.148(-e)°®°(l -b 0.37e) e<0 

For d = 7 we estimate Kc = 0.0777086(8), this is again 
closer to the series based estimate from [l9| than the MC- 
estimates from [l^, and C{Kc,oo) = 0.53(2). The limit 
specific heat is 

{ 1.75 - 2.47e°'®°(l - 0.56e -b O.Me^) ^ > o 
0.53 £ = 0 (17) 

0.064 - 0.12(-e)°'75(2+o.45e) e<0 

In both cases the uncertainty in the coefficients is in the 
last stated digit. Combining d = 5, 6, 7 and the complete 
graph case we plot them all in Figure [T^l Inside the 
scaling window, that is, with respect to k = '/N{K — 
Kc)/Kc, we can also clearly see how the specific heat for 
finite-dimensional systems approach the complete graph 
limit case. In Figure [T3] we plot C(k, L) for several linear 
orders L for d = 5,6, 7 and the complete graph. 
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FIG. 12: (Colour on-line) The limit specific heat C{e) with 
£ = {K — Kc)/Kc for d = 5 (black), d = 6 (blue) and d = 7 
(purple) and the complete graph case (red), trending down¬ 
wards at the y-axis, as do the points at e = 0. 
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FIG. 13: (Colour on-line) The specific heat C(«;) with k = 
'/N{K — Kc)/Kc- Trending downwards we see d = 5 (black, 
L > 16), d = 6 (blue, L > 8) and d = 7 (purple, L > 6) and 
the complete graph limit (red). 


DISCUSSION 

We have estimated the critical behaviour of the specific 
heat of the 5-dimensional case and derived the limit curve 
for the complete graph. The singular exponents for d = 5 
were found to differ for the high- and low-temperature 
case. To summarise, for e = {K — Kc)/Kc we estimate 
in the 5d case that for L —>■ oo the specific heat behaves 
as 

{ 2.040- 2.58e°®O(l-0.36£) e>0 

0.724 e = 0 (18) 

0.1697 - 0.231(-£)°'40(1 -b 0.26£) £ < 0 













The singular exponents are thus 0+ = 0.60 and 9~ = 
0.40 for d = 5. As d —>■ oo we expect these exponents 
to approach those of the complete graph where we find 
= 1. The exact series expansion of the limit 
specific heat for the complete graph is 

+ e>0 

C = 0.40729006421665228... e = 0 (19) 

[o £ < 0 

An open question which would be interesting to set¬ 
tle is how the left and right hand limits of the specific 
heat for the d-dimensional Ising model scales. We expect 
the limits to approach those of the mean field model as 
d —^ oo but we do not yet know how it approaches those 
values. That the d —> oo limit of the value exactly at 
should be the same as the mean field value is far from 
obvious and would also be worth further investigation. 
Similarly we would like to know the scaling with d of 
the left and right singular exponents. We expect both of 
them to approach 1, but in which way? 
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